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Coupling Constants for Scalar Glueball Decay* 
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^ ' We evaluate the partial decay widths of the lightest scalar glueball to pairs of pseudoscalar quark-antiquark 

D I states. The calculation is done in the valence (quenched) approximation on a 16^ lattice at /3 = 5.7. These 

Mh ■ predictions and values obtained earlier for the infinite volume continuum limit of the scalar glueball's mass are 

l/-^ ' in good agreement with the observed properties of /./(ITIO) and inconsistent with all other observed meson 

^SJ . resonances. 
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1. Introduction 

Whether or not glueballs have been observed 
in experiment is still generally considered unset- 
tled. Since the properties of glueballs are not ex- 
pected to be drastically different from the prop- 
erties of flavor singlet bosons including valence 
quarks and antiquarks, the identification in ex- 
periment of states virith large glueball contribu- 
tions is difficult without a reliable calculation of 
the predictions of QCD. The lattice formulation 
of QCD appears to us to give the most reliable 
method now available for determining QCD's pre- 
dictions for the masses and decay couplings of 
glueballs. 

For the infinite volume continuum limit of the 
valence (quenched) approximation to the light- 
est scalar glueball mass, we reported some time 
ago H the value 1740(71) MeV. This prediction 
was obtained using ensembles of 25000 to 30000 
gauge configurations on each of several different 
lattices. An earlier independent valence approxi- 
mation calculation of the scalar glueball mass, 
extrapolated to the continuum limit following 
Ref. §, yields a prediction of 1625(94) MeV. The 
calculation of Ref ||] uses ensembles of between 
1000 and 3000 configurations on several different 
lattices. Combining the two mass calculations 
and taking into account the correlation between 
their statistical uncertainties arising from a com- 



*Talk presented by D. Weingarten 

tPermanent address: Department of Mathematics, Trinity 

College, Dublin 2, Republic of Ireland 

tPresent address: Piazza Giovanetti 1, Novara, 28100 

Italy 



mon procedure for converting lattice quantities 
into physical units gives a scalar glueball mass 
of 1707(64) MeV. This resuh and the mass pre- 
diction with larger statistical weight are both in 
good agreement with the mass of /j(1710) and 
are strongly inconsistent with all but /o(1500) jj] 
among the other established flavor singlet reso- 
nances which could be scalars. For /o(1500) the 
disagreement is still by more than three standard 
deviations. 

The valence approximation, in effect, replaces 
the momentum and frequency dependent color 
dielectric constant arising from quark-antiquark 
vacuum polarization with its zero-momentum, 
zero- frequency limit For the long distance 
properties of hadrons, we would expect this ap- 
proximation to be fairly reliable. The infinite vol- 
ume continuum limits, for example, of the valence 
approximation to the masses of eight low lying 
hadrons composed of quarks and antiquarks dif- 
fer from experiment by amounts ranging up to 
6% A 6% error in the glueball mass would be 
100 MeV. An adaptation of an argument giving a 
negative sign for the valence approximation error 
in ^ also suggests a negative sign for the glue- 
ball mass error. Thus we would expect the scalar 
glueball in full QCD to lie between 1707(64) MeV 
and 1807(64) MeV, again favoring /j(1710) with 
/o(1500) still possible but improbable. 

The key question not answered by the mass re- 
sults, however, is whether the lightest scalar glue- 
ball has a decay width small enough for this par- 
ticle actually to be identified in experiment. It 
seems likely to us that a scalar glueball with a 
width of a few hundred MeV or less and mass 
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in the neighborhood of 1700 MeV would aheady 
have been seen in experiment. Ahernatively, if 
the scalar glueball has a width of a GeV or more, 
the prospect of ever finding this particle seems re- 
mote. A further question in the identification of 
/j(1710) as a glueball is raised by the argument 
that since glueballs are flavor singlets they should 
have the same couplings to 27ro, to 2Kl, and to 
277. This equality is violated by //(ITIO) decay 
couplings. 

In the present article we report the first lattice 
QCD calculation of the valence (quenched) ap- 
proximation to the decay couplings of the light- 
est scalar glueball to pairs of pseudoscalar quark- 
antiquark states. The couplings constants we 
obtain, combined with the mass prediction of 
1740(71) MeV, give a total two-pseudoscalar de- 
cay width of 108(29) MeV. With any reasonable 
guess concerning the scalar glueball's branching 
fraction to multibody decay modes, the result- 
ing total decay width is well below 200 MeV and 
therefore small enough for the scalar glueball to 
be identified in experiment. In fact, the predicted 
total two-pseudoscalar decay width, and individ- 
ual couplings to 27ro, to 2Kl^ and to 277 are all in 
good agreement with properties of /j(1710) and 
inconsistent with all other established flavor sin- 
glet resonances which could be scalars. The total 
two-body width of /j(1710) is 99(15) MeV [|. A 
comparison of our predicted coupling constants 
with data for /j(1710) ||^ is shown in Figure ^ 

The calculation we present uses 10500 indepen- 
dent gauge configurations on a lattice 16'^ x 24 
at (3 = 5.70. The corresponding inverse lattice 
spacing is 1.35 GeV. We believe this lattice has 
spacing sufficiently small and volume sufficiently 
large to give partial widths within 30% of their 
infinite volume continuum limits. 

In the valence approximation, according to one 
point of view, glueballs are pure field and are not 
mixed with states including valence quarks or an- 
tiquarks. The agreement between the glueball 
mass and decay couplings found in the valence 
approximation and the observed mass and decay 
couphngs of /j(1710) appears to us to be strong 
evidence that this state is largely a scalar glue- 
ball with at most some relatively smaller ampli- 
tude for configurations including valence quark- 
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Figure 1. Couplings for decay to a pair of 
pseudoscalars as a function of pseudoscalar mass 
squared. 

antiquark pairs. 

The calculations presented here were carried 
out on the GFll parallel computer ^ at IBM 
Research and took approximately two years to 
complete at a sustained computation rate of be- 
tween 6 and 7 Gfiops. A preliminary version of 
this work is discussed in Ref. [p^ . 

In the remainder of this paper we describe our 
method for determining scalar glueball decay cou- 
plings then present our numerical results. 

2. Method 

We work with a euclidean lattice gauge theory 
on a lattice x T, with the plaquette action for 
the gauge field and the Wilson action for quarks. 
We assume initially exact fiavor SU(3) symmetry 
for the quark mass matrix. Each gauge configura- 
tion is fixed to Coulomb gauge. We then define a 
collection of smeared fields. We describe smear- 
ing only for the particular choice of parameters 
actually used in the decay calculation. Let Ui{x) 
for a space direction i = 1, 2, 3, be a smeared link 
field [0 given by the average of the 9 links in direc- 
tion i from the sites of the (3 site) x (3 site) square 
oriented in the two positive space directions or- 
thogonal to i starting at site x. Let Vij{x) be the 
trace of the product of Ui{x) and Uj(x) around 
the outside of a (3 link) x (3 link) square. Let 
g{i)^ carrying zero momentum, be the sum of the 
Vij{x) for all i,j and x with time component t. 
Define ^'(a;) and ^!{x) to be quark and antiquark 
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fields, respectively, smeared by convoluting the 
local fields with a space direction gaussian with 
mean-square radius 6.0. The smeared field TTi{x) 
with flavor index i is {x)'-f^ Ai'i' (x) , where is 
a Gell-Mann flavor matrix. Define TCi {k, t) to be 
the Fourier transform of T:i{x) on the time t lat- 
tice hyperplane. 

El and are defined as the energy of a single 
pseudoscalar at rest or with momentum magni- 
tude |fc| = Stt/L, respectively. The field strength 
renormalization constant TyJ^ is defined by setting 
the large t asymptote of the vacuum expectation 
value < n\{0,t)^i{Q,Q) > to {r]^)^L^exp{-E^t). 
Define t^J similarly from a pseudoscalar field with 
momentum magnitude |fc| — 2tt/L. The glueball 
mass E^ and field strength renormalization con- 
stant 7]^ are defined by setting the large t asymp- 
tote of < git)giO)) > to be {r]3 fL^exp{-EH). 

Let the flavor singlet, rotation- 
ally invariant, two-pseudoscalar field n(ii,i2) be 
(16)"-^/^ 7ri(0, ti)7ri(0, 12), with the sum over i 
running from 1 to 8. Define the fiavor singlet field 
for two zero momentum pseudoscalars ni(ii,t2) 
as (16)^1/2 X;»'^'*(0>^i)^i(0,i2)- Let the field 
112(^1, for two pseudoscalars with opposite mo- 
menta be {2^)~'^/'^Y.^^%i{k,ti)TTi{-k,t2) where 
the sum for k is over the three positive orienta- 
tions with |fc| = 2tt/L. 

Define |1 > and |2 >, respectively, to be the 
lowest and second lowest energy flavor singlet, ro- 
tationally invariant two-pseudoscalar states, both 
normalized to 1. E^"^ is the energy of |i >. 
The amplitude rj^j^it) is then defined as L^^ < 
i\llj{t,0)\n >. At large t, ■q^J'it) approachs 
r]JJ^exp{—EJt). We expect the diagonal coeffi- 
cients rjii and ryj^ to be larger than the off- 
diagonal r]2i ^iid V12 ; respectively. Since pairs of 
pseudoscalars can interact and scatter, however, 
the off-diagonal coefficients need not be zero. 

The connected three-point functions Ti{tg,t^) 
from which we extract coupling constants are 
< g{tg)Il,{U,0) > - < g{tg) >< n,(t^,0) >. If 
the hopping constant k is chosen so that E^^ = 
E^, the lightest intermediate state which can ap- 
pear between the glueball and pseudoscalars in 
Ti{tg,tTr) is |1 >. For large tg with fixed, 
Ti{tg,tT^) will therefore be proportional to the 



coupling constant of a glueball to two pseu- 
doscalars at rest. If the k is chosen so that 
iJJ'^ = E^, however, the lightest intermediate 
state which can appear between the glueball and 
pseudoscalars in T2(ig,i7r) is still |1 >, not |2 >, 
barring the unlikely occurrence of ?yi2^(i) — 0. 
To get the coupling of a glueball to two pseu- 
doscalars with momenta of magnitude 2ttL~^ 
from T2{tg, tj^), the contribution to T2{tg, t^) from 
the |1 > intermediate state must be removed. 
We therefore define the subtracted amplitudes 
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Si{tg, t^r) = Ti{tg, t^) -———Tjltg, tjr ) , 

for (ij) of either (1,2) or (2,1). In S^{tg,U) 
the contribution from the intermediate state \ j > 
has been canceled. Although the subtraction in 
Si(tg,tT^) is irrelevant for large enough tg, we 
expect the subtracted 5i(tg,t7r) to approach its 
large tg behavior more rapidly than Ti{tg, i^) ap- 
proaches its large tg limit. 

Another state which can also appear between 
the pseudoscalars and glueball in Ti{tg,tTr) con- 
sists of a quark and an antiquark bound as a 
scalar fiavor singlet. For the lattice size, /3 and 
K used in the present calculation, this state we 
have found has a mass in lattice units above 1.25. 
The scalar glueball mass is 0.972(44). The scalar 
quark- antiquark state therefore will make only its 
appropriate virtual contribution and does not re- 
quire an additional correction. 

For tg and t.^ large with {tg + tj,) <C T, the 
three-point functions become 

■^t[tg,t^)^ S,(tg,trr). [Z) 

Here Ci = 1/V2, C2 = VS, r is (?7r2"'7Ji")/('7n ^22") 
and Ai and A2 are the glueball coupling constants 
to a pair of pseudoscalars at rest or with momenta 
of magnitude 27rL~^, respectively. As discussed 
earlier, 77^^ is defined from rj^-'[t) as large t. The 
factor Si{tg,tT^) is 



t 

ETlt- 



exp[-E9\t-tg\-E^\t\ 



U\~S,it,t^)\t-t^\]. (3) 

For t > t^r, Si{t,tTr) is the binding energy E^^ — 
2E^ and otherwise 5i{t,t.^) is 0. 
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The coupling constants in Eq. have been 
identified by fitting Si{tg, t^) with the three-point 
function arising from a simple phenomenological 
interaction lagrangian. This procedure for deter- 
mining Xi is correct to leading order in A^. A 
corresponding relation has been used for some 
time to find coupling constants among hadrons 
containing quarks and recently has produced a 
variety of results in good agreement with experi- 



ment 11 1 . The normalization of the Xi is chosen 
so that in the continuum limit they become, up 
to a factor of — i, Lorentz-invariant decay ampli- 
tudes with the standard normalization convention 
used in the section on kinematics of the Review 
of Particle Properties. 

To extract Xi from Si{tg,tTr) using Eq. we 
need the 77^^. These we obtain from propagators 
for two-pseudoscalar states. Define Ci{ti,t2) to 
be < Il{ti + 2t2,h + t2)fii{t2,0) >■ For moder- 
ately large values of ti, we then have 



aiti,t2) = Cu{t2)exp{~-Erti) 
C2^{t2)exp{-E^tl) 

c.Ah) = v7i''{t2KJ'{t2) + 



(4) 



A(t2), 



(5) 



The term Di{t2), independent of ti, arises from 
propagation across the lattice's periodic time di- 
rection boundary. Di{t2) makes a significant con- 
tribution to Ci(ti, 12) only if ^1+^2 is comparable 
to T. From Eqs. (|) and (|), 77"/ (t) and ■q^^'' can 
be found. 

One of the jobs performed by the r/^'^ in Eq. (§) 
is to correct for the interaction between the two 
pseudoscalars produced by a glueball decay. In 
the valence approximation pseudoscalars interact 
without the production and annihilation of vir- 
tual quark-antiquark pairs. The Ci{ti,t2) from 
which the vj^-^ are obtained should therefore be 
evaluated from quark propagators including only 
terms in which all initial quarks and antiquarks 
propagate through to final quarks and antiquarks. 
It can be shown that including also in the two- 
pseudoscalar propagator terms in which initial 
quarks propagate to initial antiquarks and final 
quarks propagate to final antiquark would lead 
to rflj' which correct Eq. (^) for processes missing 
from glueball decay in the valence approximation. 




Figure 2. Effective masses determined from the 
ti dependence of Ci(ti, ^2)— -Di(i2) for k = 0.1675 
and t2 — 2. 

If ti is made very large, the Ci{ti, 12) are given 
by a sum of two terms each of which is a slightly 
more complicated version of one of the exponen- 
tials in Eq. (|^) . This complication occurs because 
in the valence approximation the exchange of a 
particle composed of quarks and antiquarks be- 
tween the pseudoscalars produced in a glueball 
decay is not iterated in the same way as in full 
QCD. The term in Eq. (|) with coefficient ^^(^2) 
requires no modification if \E]'' ~ 2EJ\'^tl/2 < 1. 
Within the intervals of we use to determine the 
rjfj', these bounds are well satisfied. In any case, 
as we will discuss below, the values we obtain 
for nf/ turn out to be close to their values for 
noninteracting pseudoscalars. A consequence is 
that the corrections due to interactions between 
the decay pseudoscalars which the rjfj' contribute 
to the predicted values of Xi are small, and our 
results are fairly insensitive to details of the two- 
pseudoscalar interaction. 

3. Results 

We now turn to our numerical results. In all of 
our numerical work, fits were done by minimizing 
found from the full correlation matrix among 
the fitted data and statistical uncertainties were 
found by the bootstrap method. 

At /3 = 5.7 on a 16^ X 24 lattice, with an 
ensemble of 10500 independent configurations, 
we determined glueball and single pseudoscalar 
energies and renormalization constants following 
Rcfs. 0] and respectively. The energies in 
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Figure 3. Effective masses determined from the 
ti dependence of 6*2(^1, i2)— £'2(^2) for k = 0.1675 
and t2 = 2. 




Figure 4. Effective masses determined from the 
ti dependence 
of C2{ti,t2) - Ci2{t2)exp{-Erti) - D2{t2) for 
K = 0.1675 and <2 = 2. 
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0.48686(163) 




E9 


0.972(44) 



Table 1 

Single particle energies in lattice units, 
lattice units are shown in Table |l|. 

On a lattice of size 16'^ x 40 we then evalu- 
ated the two-pseudoscalar propagator Ci(ti,t2) 
at K = 0.1650 using 107 independent configu- 
ration, and at K — 0.1675 using 875 indepen- 
dent configurations. Fitting the ti dependence 
of Ci{ti,t2) to Eqs. (H) and (||), we determined 
and r]^J'it2) for a range of different t2- For 
each value of k we first found E^"^ by fitting 
the ti dependence of Ci{ti,t2) with fixed ^2 to 
Eq. (Q). This fit was done at large enough ti 
for the contribution of the more rapidly falling 
term C2i{t2)exp(-E^'^ti) in Eq. (|) to be neg- 
ligible. Since the coefficient C2i(t2) is signifi- 
cantly smaller than Cii(t2), it is not difficult to 
find a range of ti within which only the Cii{t2) 
term survives. For both values of n we found 
11 — ^2 < ^1 < 18 — 12 to be satisfactory. Figure |^ 
shows the fit to Ci(ti, 12) for k, — 0.1675, ^2 = 2. 
The horizontal line is the fitted value of E'^'^ . The 
points with error bars are effective energies E^f f 
given by log[ci(ii, i2)/ci(ti + 1,^2)] for Ci{ti,t2) 
defined as the propagator with constant term re- 
moved Ci{ti,t2) - Di(t2). 

We then looked for a region at large ti in 
which C2{ti,t2) gets negligible contribution from 
C22{t2)exp{-E^''ti) in Eq. (|). This region is 
more difficult to locate than the corresponding 
region for Ci(<i,<2), since the coefficient Ci2{t2) 
of the dominant large ti term in Eq. (^) is smaller 
than the coefficient 6*22(^2) of the term which falls 
more rapidly at large ti. In these fits, we there- 
fore took as input the value already determined 
for E^^ from the fit to Ci{ti,t2)- A satisfactory 
range for both k we found to be 14 — ^2 < < 
19 - t2. For the fit to (72(^1,^2) with n = 0.1675, 
t2 = 2, Figure m shows iJf in comparison to effec- 
tive energies determined from 02(^1, ^2) defined as 
C2{ti,t2) ~ D2(t2). With 6*12(^2) determined, we 
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K 


0.1650 




0.9076(49) 




Er 


1.2180(155) 


0.1675 




0.5690(44) 






0.8925(243) 




0.1650 


0.988(30) 


0.091(8) 




-0.087(8) 


1.065(13) 


0.1675 


1.050(21) 


0.107(6) 




-0.112(8) 


1.053(53) 



Table 2 

Energies in lattice units and field-strength renor- 
malizations for two-pseudoscalar states, 
then formed the subtracted propagator d2(ti,i2) 
given by C2 (i 1 , i2 ) - Ci2 {t2)exp{-Erti ) - D2 (t2) 
and looked for a plateau in effective masses at 
smaller ti to be fit to C22{t2)exp{-E^'^ti). A 
satisfactory range we found to be the 4 successive 
time values in the interval 2 < ti < 7 giving the 
fit with the smallest per degree of freedom. For 
the fit to ^2(^1,^2) with k = 0.1675, t2 = 2, Fig- 
ure ^ shows determined from 3 < ti < 6 in 
comparison to effective energies determined from 

d2{tl,t2). 

At K = 0.1650 we obtained results for < i2 < 
4, and at 0.1675 we found results for t2 of and 
2 < ^2 < 5. The values of E^'^ were statisti- 
cally consistent with being independent of t2 in 
all cases. The r]fj'{t2) were consistent with the 
asymptotic form 'i]fJ'exp{—Ejt2) in all cases for 
^2 > 2. The final values of Ef^ and ij^J' used in 
finding the were taken from combined fits with 
2 < i2 < 4 for K = 0.1650 and with 2 < ^2 < 5 
for K = 0.1675. Table | gives the final Ef'' in 
lattice units and f/fj'. For noninteracting pseu- 
doscalars fj^J^ is 1 for i = j and otherwise. Our 
data is close to these values. The final value of Ai 
is changed by less than 1 standard deviation and 
the final A2 is changed by less than 2 standard 
deviations if we ignore the determination of 77^^ 
and simply use the the noninteracting values. 

From our fitting procedure we were unable to 
obtain a satisfactory signal for the coefficient 
C2i(t2) which enters the determination of fi2i 
and therefore of Ai. Since the other -qfj^ are 
close to their values for pairs of noninteract- 
ing pseudoscalars, we took jyjj^ from the rela- 




Figure 5. Ai for tg — t.„ 
compared with a fit on tg 
6 



as a function of 
- = 0, 3 < < 7. 




Figure 6. Ai for tg — = 1 function of It^, 



compared with a fit on t„ 



1, 3 < < 7. 



tion f]2i /f]22 = —ffl-z /fill which follows from 
first order perturbation theory in the strength of 
the two-pseudoscalar interaction. This approxi- 
mation should introduce only a small error since 
simply setting 1)21 to alters the final Ai by less 
than a standard deviation. As discussed earlier, 
if Ai is determined from Eq. (H) at large enough 
ig ~ i-TTi the result is completely independent of 
V21 ■ 

From our 10500 configuration ensemble on a 
16^^ X 24 lattice, we evaluated Si and S2 for glue- 
ball decay on mass shell at k of 0.1650 and 0.1675, 
respectively. We obtained statistically significant 
results for < ig - < 2 with < < 9 for Si 
and with < < 7 for 5*2. At each point within 
this range we then determined effective A^ using 
Eq. (|^). We found Ai and A2 statistically consis- 
tent with being constant for t^^ > 3 and t^r > 2, 
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Figure 7. Ai for tg — ~ 2 as a function of tj^, 
compared with a fit on — = 2, 3 < < 7. 
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Figure 8. Ai fitted on intervals 3 < < 7 with 

a single tg — t-^ as a function of tg — t-^, compared 
with a fit on 1 < f„ - < 2. 3 < < 7. 




Figure 9. A2 for tg — 1-,^ = Q as a function of t^^, 
compared with a fit on — = 0, 2 < < 6. 




Figure 10. A2 for — = 1 as a function of It^, 
compared with a fit on f,i — t^j = 1. 2 < < 6. 




Figure 11. A2 for tg — t^ = 2 as a function of t^-, 
compared with a fit on tg — = 2, 2 < < 6. 




Figure 12. Ai fitted on intervals 3 < t^r < 7 with 
a single tg — t-j^ as a function of tg — t^, compared 
with a fit on < tg — t,r < 1, 2 < t^ < 6. 
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respectively, and all values of tg — t^^. Figures ^ 
shows effective Ai in units of the p mass as a func- 
tion of for tg — t„ — Q in comparison a fit with 
3 < < 7, — = 0. Figures || and ^ show 
corresponding data for tg — of 1 and 2, respec- 
tively. Figure ^ shows fitted values of Ai on the 
interval 3 < t,r < 7 for fixed tg — t-^ of 0, 1 or 2. 
Figures ^ shows effective A2 in units of the p mass 
as a function of tj^ for tg — t.^ = in comparison 
a fit with 2 < t^ < 6, tg — tT^ = 0. Figures |l^ and 
show corresponding data for tg — of 1 and 
2, respectively. Figure |l^ shows fitted values of 
Ai on the interval 2 < t.^ < 6 for fixed tg — of 
0, 1 or 2. 

To extract final values of Ai, we tried fits to all 
rectangular intervals of data including at least 4 
values of and at least 2 values of tg — t-^. For 
each Ai we chose the fit giving the lowest value 
of per degree of freedom. The window de- 
termined in this way for Ai is 3 < ^ 7 with 
1 < tg — t„ < 2, and for A2 is 2 < < 6 with 
< tg — t^r < 1. The horizontal lines in Fig- 
ures II and |l2| shows the final value of Ai and A2 
respectively. Over the full collection of windows 
we examined, the fitted results varied from our 
final results by at most 1 standard deviation. We 
believe our best fits provide reasonable estimates 
of the asymptotic coefficients in Eq. (||). 

So far our discussion has been restricted to 
QCD with u, d and s quark masses degenerate. 
An expansion to first order in the quark mass ma- 
trix taken around some relatively heavy SU(3) 
symmetric point gives glueball decay couplings 
for tt's, K's and the 77 which are a common lin- 
ear function of each meson's average quark mass. 
Since meson masses squared are also nearly a lin- 
ear function of average quark mass, the decay 
couplings are a linear function of meson masses 
squared. Thus from a linear fit to our predic- 
tions for decay couplings as a function of pseu- 
doscalar mass squared at unphysical degenerate 
values of quark masses we can extrapolate de- 
cay couplings for physical nondegenerate values 
of quark masses. From this linear fit a predic- 
tion can also be made for the decay coupling of 
the scalar glueball to t] + ri' , if we ignore the con- 
tribution to the decay from the process in which 
the rj quark and antiquark are connected to each 



other by one propagator and the 77' quark and an- 
tiquark are connected to each other by a second 
propagator. 

Figure |l| shows predicted coupling constants as 
a function of predicted meson mass squared along 
with linear extrapolations of the predicted values 
to the physical tt, K and rj masses, in compari- 
son to observed decay couplings |^ for decays of 
/j(1710) to pairs of tt's, K's and ry's. Masses and 
decay constants arc shown in units of the p mass. 
Our predicted width for the scalar glueball decay 
to 1] + rj' is 6(3) MeV. For the ratio A,p,'/A,,^ we 
get 0.52(13). We predict a total width for glue- 
ball decay to pseudoscalar pairs of 108(29) MeV, 
in comparison to 99(15) MeV for /j(1710). 

One of us (D.W.) is grateful to S. Lindenbaum, 
R. Longacre, S. Sharpe and the participants in 
Gluonium 95 for valuable conversations. 
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